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ABSTRACT 

An analytical method has been developed 

to determine gear tooth surface coordinates of 
face-milled spiral bevel gears. The method com- 
bines the basic gear design parameters with the 
kinematical aspects for spiral bevel gear manu- 
facturing. A computer program was developed to 
calculate the surface coordinates. From this data 
a three-dimensional model for finite element ana- 
lysis can be determined. Development of the 
modeling method and an example case are presented. 

INTRODUCTION 

Spiral bevel gears are currently used in all 
helicopter power transmission systems. This type 
of gear is required to turn the corner from a 
horizontal engine to the vertical rotor shaft. 
These gears carry large loads and operate at high 
rotational speeds. Recent research has focused on 
understanding many aspects of spiral bevel gear 
operation, including gear geometry [1-121 , gear 
dynamics [13-15] , lubrication [16] , stress analy- 
sis and measurement [17-21] , misalignment [22,22] , 
and coordinate measurements [24,25] , as well as 
other areas. 

Research in gear geometry has concentrated on 
understanding the meshing action of spiral bevel 
gears [8-11 ] . This meshing action often results 
in much vibration and noise due to an inherent 
lack of conjugation. Vibration studies [26 ] have 
shown that in the frequency spectrum of an entire 
helicopter transmission, the highest response can 
be that from the spiral bevel gear mesh. There- 
fore if noise reduction techniques are to be im- 
plemented effectively, the meshing action of 
spiral bevel gears must be understood. 

Also, investigators [18,19] have found that 
typical design stress indices for spiral bevel 
gears can be significantly different from those 
measured experimentally. In addition to making 
the design process one of trial and error (forcing 
one to rely on past experience), this inconsis- 
tency makes extrapolating over a wide range of 
sizes difficult, and an overly conservative design 
can result. 

The objective of the research reported herein 
was to develop a method for calculating spiral 
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bevel gear-tooth surface coordinates and a three- 
dimensional model for finite element analysis. 
Accomplishment of this task requires a basic un- 
derstanding of the gear manufacturing process, 
which is described herein by use of differential 
geometry techniques [ 1 ] . Both the manufacturing 
machine settings and the basic gear design data 
were used in a numerical analysis procedure that 
yielded the tooth surface coordinates. After the 
tooth surfaces (drive and coast sides) were de- 
scribed, a three-dimensional model for the tooth 
was assembled. A computer program was developed 
to automate the calculation of the tooth surface 
coordinates, and hence, the coordinate for the 
gear-tooth three-dimensional finite element model. 
The basic development of the analytical model is 
explained, and an example of the finite element 
method is presented. 

DETERMINATION OF TOOTH SURFACE COORDINATES 

The spiral gear machining process described 
in this paper is that of the face-milled type. 
Spiral bevel gears manufactured in this way are 
used extensively in aerospace power transmissions 
(i.e., helicopter main/tail rotor transmissions) 
to transmit power between horizontal gas turbine 
engines and the vertical rotor shaft. Because 
spiral bevel gears can accommodate various shaft 
orientations, they allow greater freedom for over- 
all aircraft layout. 

In the following sections the method of de- 
termining gear-tooth surface coordinates will be 
described. The manufacturing process must first 
be understood and then analytically described. 
Equations must be developed that relate machine 
and workpiece motions and settings with the basic 
gear design data. The simultaneous solution of 
these equations must be done numerically since no 
closed-form solution exists. A description of 
this procedure follows. 

Gear Manufacture 

Spiral bevel gears are manufactured on a 
machine like the one shown in Fig. 1. This 
machine cuts away the material between the concave 
and convex tooth surfaces of adjacent teeth simul- 
taneously. The machining process is better illus- 
trated in Fig. 2. The head cutter (holding the 
cutting blades or the grinding wheel) rotates 
about its own axis at the proper cutting speed, 
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Figure 2.— Orientation of workpiece to generation machinery. 


about its own axis at the proper cutting speed, 
independent of the cradle or workpiece rotation. 
The head cutter is connected to the cradle through 
an eccentric that allows adjustment of the axial 
distance between the cutter center and cradle 
(machine) center, and adjustment of the angular 
position between the two axes to provide the de- 
sired mean spiral angle. The cradle and workpiece 
are connected through a system of gears and 
shafts, which controls the ratio of rotational 
motion between the two (ratio of roll). For 
cutting, the ratio is constant, but for grinding, 
it is a variable. 

Computer numerical controlled (CNC) versions 
of the cutting and grinding manufacturing pro- 
cesses are currently being developed. The basic 
kinematics, however, are still maintained for the 
generation process; this is accomplished by the 
CNC machinery duplicating the generating motion 
through point-to-point control of the machining 
surface and location of the workpiece. 

Coordinate Transformations 

The surface of a generated gear is an 
envelope to the family of surfaces of the head 
cutter. In simple terms this means that the 
points on the generated tooth surface are points 
of tangency to the cutter surface during manu- 
facture. The conditions necessary for envelope 
existence are given kinematically by the equation 
of meshing. This equation can be stated as fol- 
lows: the normal of the generating surface must 

be perpendicular to the relative velocity between 
the cutter and the gear-tooth surface at the point 
in question [ 1 ] . 

The coordinate transformation procedure that 
will now be described is required to locate any 
point from the head cutter into a coordinate 
system rigidly attached to the gear being manu- 
factured. Homogeneous coordinates are used to 


allow rotation and translation of vectors simply 
by multiplying the matrix t rans f ormat ions . The 
method used for the coordinate transformation can 
be found in Refs. 1, 5, 8 to 11, and 27. 

To begin the process we start with the head 
cutter (Fig. (3)). In this report it is assumed 
that the cutters are s t raight- s ided . The para- 
meters u and 6 determine the location of a 
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current point in the cutter coordinate system S c . 
Angles and ^ ob are the blades that cut the 

convex and concave sides of the gear tooth respec- 
tively. Thus 7 c is given by the following 
equation : 

r cot ^ - u cos ^ 
u sin ^ sin 9 ( 1 ) 

c u sin i/) cos 9 
1 

Once we have 7. we then transform from one coor- 
dinate system to the next for the coordinate sys- 
tems as shown in Figs. (4) and (5) [27] . 

Using matrix transformations we can determine 
the coordinates in S w of a point on the generat- 
ing surface by: 


Here [M i .] describes the required homogeneous 
coordinate transformation from system " j " to 
system "i". Therefore Eq. (2) describes the 
location of a point in the gear fixed coordinate 
system based on machine settings L n , E n , q, s, r, 
, and gear design information ji and 6. At 
this point the machine settings and the gear 
design values are known. Parameters u. 9, and 
<p c are the unknown variables that are solved for 
numerically . 




Figure 5.— Top and front view of right hand gear surface generation 
(Oc = 0 shown here). 


Tooth Surface Coordinate Solution Procedure 

In order to solve for the coordinates of a 
spiral bevel gear-tooth surface, the following 
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Figure 4.— Orientation of cutter, cradle, and fixed coordinate 
systems. S c . S s . and S m respectively. 


items must be used simultaneously: the trans- 

formation process, the equation of meshing, and 
the basic gear design information. The trans- 
formation process described previously is used to 
determine the location of a point on the head cut- 
ter in coordinate system S w . Since there are 
three unknown quantities (u, 9, and <p c ) , three 
equations relating them must be developed. 

Values for u, 9 , and <p c are used to satisfy 
the equation of meshing given by Refs. 1 and 9. 

n • V = O ( 3 ) 

where n is the normal vector to the cutter and 
workpiece surfaces at the specified location of 
interest, and V is the relative velocity between 
the cutter and workpiece surfaces at the specified 
location . 

Gear design information is then used to 
establish an allowable range of values of the 
radial (?) and axial (z) positions that are known 
to exist on the gear being generated. This is 
shown in Fig. 6 . 

First the equation of meshing must be satis- 
fied. This is given as [ 9 1 : 

(u - rcot cos COS 7 sin ; q 1 <p c j 

+ S - sin7jcos ^sin# * cos7sin^sin|q - 0 c jj 

± E m |cos 7 sin^ + sin 7 cos^cos |0 * q ± ^ c jJ 
- L m | sin 7 cos ^sin |0 * q ± <p c jj - 0 
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Figure 6.— Orientation of gear to be generated with assumed positions 
r and z. 


In summary the procedure is as follows. 

Known locations r and z on the active profile 
are used along with the equation of meshing to 
determine unknown parameters u, d t and (p . The 
parameter values, machine tool settings, and gear 
design values are used in the coordinate trans- 
formation shown earlier to find the radial and 
axial positions in the gear fixed coordinate 
system . This procedure results in the solu- 
tion of three simultaneous nonlinear algebraic 
equations that are solved numerically. 

APPLICATION OF SOLUTION TECHNIQUE 

An application of the techniques previously 
discussed will now be presented. The component to 
be modeled was from the NASA Lewis Spiral Bevel • 
Gear Test Facility. A photograph of the spiral 
bevel gear mesh is shown in Fig. 8, and the design 
data for the pinion member are shown in Table 1. 


The upper and lower signs preceding the 
above terms pertain to left and right hand gears 
respectively . 

The axial position must match the value found 
from transforming the cutter coordinates S c to 
workpiece coordinates S w - This is satisfied by 
the following (Fig. 6): 


z w - z = 0 ^ 

Finally the radial location from the work 
axis of rotation must be satisfied. This is 
accomplished by using the magnitude of the loca- 
tion in question in the x w - y w plane (Fig. 6): 



Now a system of three equations (Eqs. (3) to 
(5)) is solved simultaneously for the three para- 
meters u, 6, and 0 c , for a given gear design with 
a set of machine tool settings. These are non- 
linear algebraic equations that can be solved 
numerically with commercially available mathe- 
matical subroutines. These equations are then 
solved simultaneously for each location of in- 
terest along the tooth flank, as shown in Fig. 7. 
From the surface grids, the active profile 
(working depth) occupied by a single tooth is 
defined . 



Figure 7.— Calculation points (100 each side) for concave and 
convex sides of tooth surface. 


Figure 8— Photograph of NASA spiral bevel 
gear test rig components. 

Surface Coordinate Calculation 

Using Figs. 7 and 9 as references, we will 
describe the calculation procedure for surface 
coordinates. First, the concave side of the tooth 
is completely defined before moving to the convex 
side. These points are calculated by starting at 
the toe end and at the lowest point of active 
profile. Nine steps of equal distance are used 
from the beginning of the active profile to the 
face angle (addendum) of the gear tooth, and then 
back to the next axial position (Fig. 9). The 
procedure is repeated until the concave side is 
completely described. Then the same procedure is 
followed for the convex side. 

Example Model and Results 

From the one-tooth model described earlier 
the analysis techniques can be demonstrated. The 
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TABLE I. 


EXAMPLE CASE FOR SURFACE COORDINATE GENERATION 


Gear 

design data 


Number of teeth (pinion, gear) 



12, 36 

Dedendum angle, <5, deg . . . . 



1.0 

Addendum angle, deg 



3.883 

Pitch angle, p, deg 



18.433 

Shaft angle, deg 



90.0 

Mean spiral angle, deg . . . . 



35.0 

Face width, mm, (in.) . . . . 



. . . . 25.4 (1.0) 

Mean cone distance, mm (in.) 



. . . 81.05 (3.191) 

Inside radius of gear blank, mm (in.) 


. . . 15.3 (0.6094) 

Top land thickness, mm (in.) 



. . . 2.032 (0.080) 

Clearance mm (in.) 



. . . 0.762 (0.030) 

Generation 

machine 

settings 



Concave 

Convex 

Radius of cutter, r, mm (in.) 

75.222 

(2.9615) 

78.1329 (3.0761) 

Blade angle, i// , deg 

161.358 


24.932 

Vector sum, L , mm, in. 

1.0363 

(0.0408) 

-1.4249 (-0.0561) 

Machine offset, E a , mm (in.) 

3.9802 

(0.1567) 

-4.4856 (-0.1766) 

Cradle to cutter distance. 




s , mm ( in . ) 

74.839 

(2.9646) 

71.247 (2.3050) 

Cradle angle, q, deg 

64.01 


53.82 


v- Surface coordinate 



1724 MPa (250 ksi) 
pressure load 
normal to gear 



Figure 10— Boundary conditions for the constant fillet/root radius model for the example application. 
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model shown in Fig. 10 is that for a constant 
fillet and root radius (also called full fillet) 
model. The fillet and root radius on the convex 
side has been added along with the tooth section 
(without the tooth) to make the model symmetric 
about the tooth centerline. Figure 10 shows a 
hidden line plot of the finite element mesh with 
eight-noded isoperimetric three-dimensional solid 
continuum elements. This model has 765 elements 
and 1120 nodes. The boundary conditions are also 


shown in Fig. 10. A 1724-MPa (250 ksi) constant 
pressure load was applied normal to the tooth 
surface of nine elements, and the two edge sur- 
faces of the gear rim had all degrees of freedom 
constrained . 

The results were calculated by MSC/NASTRAN 
and were subsequently displayed by PATRAN. Fig- 
ure 11 shows the major principle stresses for the 
boundary conditions shown in Fig. 10. 
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-34 000 

A 

-234 
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14 000 

E 


26 000 

F 
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H 


62 000 

1 
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J 
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86 000 

K 


98 000 

L 


110 000 

M 

759 


Figure 1 1 . — Major principle stress for the boundary conditions specified in figure 1 0. 


SUMMARY OF RESULTS 

A method has been presented that uses dif- 
ferential geometry techniques to calculate the 
surface coordinates of face-milled spiral bevel 
gear teeth. The coordinates must be solved for 
numerically by a simultaneous solution of non- 
linear algebraic equations. These equations 
relate the kinematics of manufacture to the gear 
design parameters. Coordinates for a grid of 
points are determined for both the concave and 
convex sides of the gear tooth. These coordinates 
are then combined to form the enclosed surface of 
one gear tooth. A computer program, was developed 
to solve for the gear-tooth surface coordinates 
and provide input to a three-dimensional geometric 
modeling program (i.e., PATRAN). This enables an 
analysis by the finite element method. An example 
of the technique was presented. 
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